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Abstract 

The "entanglement of formation" of a mixed state of a bipartite quantum 
system can be defined in terms of the number of pure singlets needed to 
create the state with no further transfer of quantum information. We find an 
exact formula for the entanglement of formation for all mixed states of two 
qubits having no more than two non-zero eigenvalues, and we report evidence 
suggesting that the formula is valid for all states of this system. 
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Entanglement is the potential of quantum states to exhibit correlations 
that cannot be accounted for classically. For decades, entanglement has been 
the focus of much work in the foundations of quantum mechanics, being asso- 
ciated particularly with quantum nonseparability and the violation of Bell's 
inequalities In recent years, however, it has begun to be viewed also 
as a potentially useful resource. The predicted capabilities of a quantum 
computer, for example, rely crucially on entanglement @, and a proposed 
quantum cryptographic scheme achieves security by converting shared en- 
tanglement into a shared secret key || . For both theoretical and potentially 
practical reasons, it has become interesting to quantify entanglement, just as 
we quantify other resources such as energy and information. In this letter we 
adopt a recently proposed quantitative definition of entanglement and derive 
an explicit formula for the entanglement of a large class of states of a pair of 
binary quantum systems (qubits). 

The simplest kind of entangled system is a pair of qubits in a pure 
but non-factorizable state. A pair of spin-| particles in the singlet state 
TTfd 1 1) — I IT)) is perhaps the most familiar example, but one can also 
consider more general states such as a\ fl) + /3\ if), which may be less en- 
tangled. For any bipartite system in a pure state, Bennett et al. Q have 
shown that it is reasonable to define the entanglement of the system as the 
von Neumann entropy of either of its two parts. That is, if \ip) is the state of 
the whole system, the entanglement can be defined as E(ip) = — Trplog 2 p, 
where p is the partial trace of over either of the two subsystems. (It 

does not matter which subsystem one traces over; the result is the same ei- 
ther way.) What Bennett et al. showed specifically is the following. Consider 
n pairs, each in the state Let an observer Alice hold one member of each 
pair and let Bob, whom we imagine to be spatially separated from Alice, hold 
the other. Then if has E "ebits" of entanglement according to the above 
definition, the n pairs can be reversibly converted by purely local operations 
and classical communication into m pairs of qubits in the singlet state, where 
m/n approaches E for large n and the fidelity of the conversion approaches 
100%. This interconvertibility is strong justification for the above definition 
of E and characterizes it uniquely. 

It is somewhat harder to define the entanglement of mixed states [0], 
though again one can use the singlet as the basic unit of entanglement and 
relate the given mixed state to singlets. The new feature in the case of 
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mixed states is that the number of singlets required to create the state is not 
necessarily the same as the number of singlets one can extract from the state 
JJ. In this paper we focus on the former quantity, which leads to the concept 
of "entanglement of formation." In mathematical terms the entanglement of 
formation is defined as follows 0. Given a mixed state p of two quantum 
systems A and B, consider all possible ways of expressing p as an ensemble 
of pure states. That is, we consider states and associated probabilities 
Pi such that 

P = ^2Pi\^i)(^i\- (1) 

i 

The entanglement of formation of p, E(p), is defined as the minimum, over 
all such ensembles, of the average entanglement of the pure states making 
up the ensemble: 

£ = min Y^PiE^i). (2) 

i 

Entanglement of formation has the satisfying property that it is zero if and 
only if the state in question can be expressed as a mixture of product states. 
For ease of expression, we will refer to the entanglement of formation simply 
as "entanglement." 

Peres || and Horodecki et al. || have found elegant characterizations 
of states with zero and non-zero E, and Bennett et al. have determined 
the value of E for mixtures of Bell states. (These are a particular set of 
orthogonal, completely entangled states of two qubits; we will refer to other 
sets of such states as "generalized Bell states.") But the value of E for most 
states, even of two qubits, is not known, and in fact it has not been evident 
that one can even express E in closed form as a function of the density 
matrix. The exact formula we derive in this letter is proved for all density 
matrices of two qubits having only two non-zero eigenvalues, but it appears 
likely that it applies to all states of this system. 

Our starting point is a curious and useful fact about the pure states of 
a pair of qubits. For such a system, we define a "magic basis" consisting of 
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the following four states (they are the Bell states with particular phases) [0]: 

|ei) = f(|n> + iU» 

|e 2 > = fz(|TT>-|U» ^ 
|e 3 > = fi([U) + UT» [) 
|e4> = |(|U)-UT» 

where we have used spin-| notation for definiteness. When a pure state \ip) 
is written in this particular basis, as \ip) = J2i a i\ e i), hs entanglement can 
be expressed in a remarkably simple way in terms of the components af 
Define the function 

£(x) = F(- + -\fl - x 2 ) for < x < 1, (4) 

where H is the binary entropy function H(x) — — [x log 2 x+(l— x) log 2 (l— x)]. 
Then the entanglement of is 

E{i/>) = £{C{i/>)), (5) 

where C is defined by 

c(V0 = l£«il- (6) 

i 

The quantity C, like -E for this system, ranges from zero to one, and it 
is monotonically related to E, so that C is a kind of measure of entangle- 
ment in its own right. It is sufficiently useful that we give it its own name: 
concurrence. As we look for a pure-state ensemble with minimum average 
entanglement for a given mixed state, our plan will be to look for a set of 
states that all have the same entanglement, which is to say that they all have 
the same concurrence. 

Two other facts about the magic basis are worth highlighting, (i) The 
set of states whose density matrices are real when expressed in the magic 
basis is the same as the set of mixtures of generalized Bell states (Horodecki 
et al. have called such mixtures "T states" [0). (ii) The set of unitary 
transformations that are real when expressed in the magic basis (or real 
except for an overall phase factor) is the same as the set of transformations 
that act independently on the two qubits. 
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It happens that our formula for E is conveniently expressed in terms of 
a matrix R, which is a function of p defined by the equation 



R(p) = ^Vpp*Vp- (7) 

Here p* is the complex conjugate of p when it is expressed in the magic basis; 
that is, p* = \ e i) ( e j\p\ e i) ( e j\- To get some sense of the meaning of R, 
note that Tr R, ranging from to 1, is a measure of the "degree of equality" 



11|] between p and p*, which in turn measures how nearly p approximates a 
mixture of generalized Bell states. Note also that the eigenvalues of R are 
invariant under local unitary transformations of the separate qubits, a fact 
that makes these eigenvalues particularly eligible to be part of a formula for 
entanglement, since entanglement must also be invariant under such trans- 
formations. We now state our main result. 

Theorem. Let p be any density matrix of two qubits having no more than 
two non-zero eigenvalues. Let X max be the largest eigenvalue of R(p). Then 
the entanglement of formation of p is given by 

E(p)=S(c); c = max(0,2X max -TiR). (8) 

(The quantity c can thus be called the concurrence of the density matrix p. 
If p is pure, then c reduces to the concurrence defined in Eq. @.) 

Proof. Let \vi) and |i; 2 ) be the two eigenvectors of p corresponding to 
its two non-zero eigenvalues. Define the 2x2 matrix r such that = 
Vi ■ Vj, where the dot product is taken in the magic basis with no complex 
conjugation: Vi-Vj = J2k ( e k\ Vi) (ek\ Vj). Consider an arbitrary pure state \ip) 
that can be written in the form = a \v\) + b \v2). If \ip) is expressed as a 
4- vector in the magic basis, we can rewrite Equation (||) as C(ip) = \ip ■ ip\, 
and 

C 2 (iP) = ■ ■ $)* = Tr [s*tst*\, (9) 

where s=^^^(a b ) is the density matrix of \ip) in the (i^i^-basis. 

Let us define the function 

f(u) = Tr [uj*tujt*\ (10) 

for any density matrix ui expressed in the (t^i^-basis. From Eq. @, f(ui) = 
C 2 (uj) if lo represents a pure state. Now, to is a 2 x 2 density matrix, and 
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as such can be written as a real linear combination of Pauli matrices: 00 = 
|(7 + r • a) where Tj = Tr [cjuj]. Substituting this form into Eq. ([T0|) gives 
us an expression 

/H = ^Tr [r*r] + £ r,L, + £ r,r,M^ (11) 

with 

% = ~Tr h-rV] (12) 

and 

M ij = ±Tr[otT<Trf. (13) 

Thus / is defined on the surface and interior of a unit sphere in three dimen- 
sions, the domain of r. 

M is a real, symmetric matrix with eigenvalues ±|| detr| and jTr [t*t], 
and L is the eigenvector of M corresponding to the eigenvalue ^Tr [r*r]. 
Since M has two positive eigenvalues and one negative eigenvalue, f{oj) is 
convex along two directions and concave along a third. For the purpose of 
this proof, we would like to have a function that is equal to f(s) for pure 
states s, but convex in all directions. With this in mind we define 

sM = /H + i|detT|(|rf-l), (14) 

which is identical to / for pure states (|rj = 1). The extra term added to 
f(u>) in effect adds a constant to / and a multiple of the identity matrix to 
M. If we define a matrix 

N = M + i|detr|J (15) 

and a constant 

K = -Tr [t*t] + ^ I det r| (16) 

then we can write 

//U) A' • • X'-'-A,- (17) 

3 h3 
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The added term in Eq. ([15|) makes all the eigenvalues of N non-negative, 
one of them being zero. Thus g is a convex function. Since L is an eigen- 
vector of N associated with a positive eigenvalue, and is orthogonal to the 
eigenvector with zero eigenvalue, the function g is constant along the latter 
direction. We can imagine the function g (suppressing one dimension) as a 
sheet of paper curved upward into a parabolic shape; it achieves its minimum 
value along a straight line. Moreover, one can show by direct calculation that 
the minimum value of g is zero. In the sphere where r takes its values, it is 
interesting to imagine the surfaces along which g is constant. These surfaces 
appear as co-axial cylinders with elliptical cross-sections. (The common axis 
of the cylinders is the line along which g = 0.) The mixed state p that we 
are considering lies on one of these cylinders and can be decomposed into 
two pure states lying on the same cylinder, that is, having the same value 
of g. The next two paragraphs show that no other decomposition of p has a 
smaller average entanglement than this one. 

Any decomposition of p into pure states can be viewed as a collection of 
weighted points on the surface of the sphere, such that the "center of mass" 
of these points is the point representing p. The average entanglement of such 



an ensemble is the average of £(Jg(s)) over the ensemble, since £(Jg(s)) 



is equal to entanglement for pure states s. If we can show that £(^g(uj)), 
regarded as a function of uo, is convex over the interior of the sphere, then 



it will follow that this average cannot be less than £(yg(p)). But we have 
just seen that p can be decomposed into two pure states s for which g(s) is 
the same as g{p), so this will prove that the entanglement of p is equal to 

In fact it is not hard to prove the desired convexity. The function g{uj) 
is parabolic with minimum value zero. Its square root is therefore a kind of 
cone and is also convex. The function £(x) is a convex and monotonically 
increasing function of x. It follows, then, from the transitive property of 



convex functions [12| that £(y g(u)) is a convex function of lu. 

We have thus found the entanglement of p and need only express it in a 
simpler form. Replacing u with p in Eq. (|T0|) and using the fact that p is 
diagonal in the (v i, t>2)-basis, we obtain 

f(p) = Tr (R 2 ) = Xj + \l (18) 
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where Ai and A2 are the non-zero eigenvalues of R (Eq. ([?[)). Similarly, one 
finds that for the other term in Eq. (14), 



-|detr| (|f| 2 - l) = -2AiA 2 , (19) 

so that g(p) = \\ + \% — 2AiA 2 . Taking the square root, we arrive at the 
result 

E(p)=S(c); C =|A 1 -A 2 |. (20) 

The expression (P0|) is equivalent to Eq. (R|) for the case of two non-zero 
eigenvalues. This completes the proof of the theorem. 

Although we have proved our result only for density matrices with just 
two non-zero eigenvalues, we can report three pieces of evidence suggesting 
that the formula fl8|) may hold quite generally for a system of two qubits. 

1. For a mixture of Bell states, mixed with probabilities pi, . . . , p$, Ben- 
nett et al. [7[ have shown that the entanglement is equal to S(c), with 
c given by max (0, 2p max — 1). But in this case R is equal to p, so that 
our expression is equal to theirs. Thus our formula applies also to this 
class of density matrices, most of which are not covered by the above 
theorem. 

2. Peres || and Horodecki et al. || have provided a test, based on partial 
transposition, for determining whether a given state of two qubits has 
zero or non-zero E. We have applied both the Peres-Horodecki test 
and our own formula to several thousand randomly chosen density ma- 
trices and have found agreement between them in every case. That is, 
Eq. (|3|) gave E > if and only if the Peres-Horodecki test indicated 
the presence of entanglement, which happened in roughly one-third of 
the cases. 

3. For each of twenty-five randomly chosen density matrices with non- 
zero entanglement, we have explored the space of all decompositions of 
the density matrix into pure states, limiting ourselves to ensembles of 
four states. (The example of Bell mixtures [0] suggests that four-state 
ensembles may be sufficient.) In each case, the result of numerically 
minimizing the average entanglement of the ensemble agrees with the 
result predicted by our formula. 
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If the formula turns out to be correct for all states, it will considerably 
simplify studies of entanglement. Questions such as whether the "distill- 
able entanglement" is equal to the entanglement of formation || [/J, that is, 
whether one can extract as much entanglement as one puts into the state, 
will presumably be easier to answer if there is an explicit formula for the 
latter quantity. It is also conceivable that our result can be generalized to 
systems with larger state-spaces, such as an entangled pair of n-level atoms, 
though it is not clear whether there is any structure in such spaces that 
would play quite the same role that the magic basis plays in the two-qubit 
case. In imagining possible generalizations, it is interesting to note that the 
form of R has much in common with the "mixed-state fidelity" Jll| of Bures, 
Uhlmann, and Jozsa, which is in no way special to two-qubit systems. 

We would like to thank Charles Bennett, David DiVincenzo and John 
Smolin for helpful and stimulating discussions. 
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